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^ ■ Abstract 

^ \ The simulations of the hght scalar mesons on the lattice are presented at the introductory 

level. The methods for determining the scalar meson masses are described. The problems 
^ ■ related to some of these methods are presented and their solutions discussed. 
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O 1 Introduction 

The observed spectrum of the light scalar resonances below 2 GeV is shown in Fig. ^ 
The existence of flavor singlet a and strange iso-doublet k are still very controversial p. 
. Irrespective of their existence, it is difficult to describe all the observed resonances by one 
or two SU (3) flavor nontes of qq states: 

^1 • If cr and K, do not exist, than J<'o(1430) has to be strange partner of ao(980), but the 

mass difference appears to big. Also there are to many states to be described by one 
^ ■ nonet. 

• If (7 and K exist, then all these states could represent two qq nonets and one glueball, 
^ ■ where the largest glueball component is commonly attributed to /o(1500). However, 

■ most of the models and lattice simulations have difficulties in relating the observed 

properties of states below 1 GeV to the qq states. 

This situation is in contrast to the spectrum of light pseudoscalar, vector and axial- vector 
resonances, where qq assignment works well. It raises a question whether the scalar reso- 
nances below 1 GeV are conventional qq states or perhaps exotic states such as tetraquarks 

m 

This issues could be settled if the mass of the lightest qq states could be reliably deter- 
mined on the lattice and identified with the observed resonances. In lattice QCD, the hadron 
masses are conventionally extracted from the correlation functions that are computed on the 
discretized space-time. 



^Talk presented at Mini-workshop Exciting Hadrons at Bled, Slovenia, July 2005. 
^Electronic address: sasa.prelovsek@ijs.si 
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Figure 1: The spectrum of observed light scalar resonances below 2 GeV The existence of a 
and K are still very controversial experimentally. 

In the next section we present how the scalar correlator is calculated on the lattice. The 
relation between the scalar correlator and the scalar meson mass is derived in Section 3. A 
result for the mass of / = 1 scalar meson is presented in Section 4. In Section 5 we point 
out the problems which arise due to the unphysical approximations that are often used in 
the lattice simulations and we discuss the proposed solutions. We close with Conclusions. 

This article follows the introductory spirit of the talk given at the Workshop Exciting 
hadrons^ and many technical details are omitted. 

2 Calculation of the scalar correlator 

Let us consider the correlation function for a flavor non-singlet scalar meson qiq2 first. In 
a lattice simulation it is calculated using the Feynman functional integral on a discretized 
space-time of finite volume and finite lattice spacing. The correlation function represents a 
creation of a pair qiq2 with = O"*" at time zero and annihilation of the same pair at some 
later Euclidean time t 



where both quarks are created (annihilated) at the same spatial point for definiteness here'^. 
Wick contraction relates this to the product of two quark propagators shown by the connected 
diagram in Fig. 
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The quark propagator in the gluon field G and Euclidean space-time |3] 





•^Different shapes of creation and annihilation operators in spatial direction can be used. 
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is the inverse of the discretized Dirac operator ^£;+mj , which is a matrix in coordinate space 
and depends on the gluon field G ^. The inversion of a large Dirac matrix is numerically 
costly, but the calculation of correlator Q is feasible since it depends only on two propagators 
from a certain point (0,0) to all points (a:, t). Both of these are obtained by solving the 
equation {]pE + ^i)y' = V for a single^ source vector V which is non-zero only at (0, 0). The 
expectation value over the gluon fields in Q is computed based on the Feynman functional 
integral 

^' JVGjVqJVq e-^QCD J VG H.det [Pe + rrii] e^^o ' ^ ^ 

A finite ensemble of gluon field configurations is generated in the lattice simulations. Each 
configuration is generated with a probability njdet[^£; + mi\ e~^^ for a given discretized 
gauge action Sg and Dirac operator pE- The functional integral (jlj) is calculated as a sum 
over the ensemble 

1 ^ 
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Figure 2: The disconnected (a) and the connected (b) Feynman diagrams that need to be evaluated 
to compute the correlator. The disconnected diagram is present only for the flavor singlet meson. 

The correlator for the flavor singlet scalar meson qq 

C{t) = 5^(0|g(f,t)g(f,t) g(0,0)g(0,0)|0) (6) 

X 

requires also the calculation of the disconnected diagram in Fig. 

Tr,,cPropg 0^0,0 X! Tr,,eProp^^,t^^,t ) (7) 

in addition to connected one. The propagator Prop^j^^j in principle requires the solution 
oi {IpE ^ 'mi)V' = V for source vector V at any point. Such a number of inversions is 
normally prohibitively large and one is forced to use approximate methods for evaluating 
the disconnected part (|7j) of the singlet correlator. A calculation of the correlator for singlet 
meson in therefore much more demanding than for non-singlet meson. 



'^P = 7^(0^ + |AaG°) in continuum Minkowski space-time. 

■''In fact (Pe + mi)V' = V has to be solved for every spin and color of the source vector V. 
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3 Relation between correlator and meson mass 



In this Section we derive the relation between the scalar correlator and the scalar meson 
mass. The state ^(0)g(0)|0) that is created at time zero is not a scalar meson 15"), but it is a 
superposition of the scalar meson and all the other eigenstates of Hamiltonian \n) with the 
same quantum numbers = 0"^ and I'^ as \qq) 

m = j: c„i„) = c,is> + CIS.) + 5: c. ^ 7'" \ + ... UcM -"';7) . (8) 

^-^ hadron st. / ^ \ smglet / 

Here IS*) and 15*) are ground and excited scalar mesons, while the third term represents the 
sum over multi-hadron states. The eigenstate In) evolves as e^P"X-^"t Euclidean space- 
time, so the scalar correlators (0) and © evolve as 

C{t) =J2{mt)q{x,t) g(0,0)g(0,0)) $^(gg|n)e^^^"^^-^"*(n|gg) =^|(gg|n)|2 e" 



■E„t I 



\p=0 



V-^* + Ic^l^e-™-* + V |q|V^^^^""* + ... +\co\' ) . (9) 

^-^ \ smglet ' 



If \S) is the lightest state among \n), than C{t) oc e~"^^* at large t and ms and can be 
extracted simply by fitting the lattice correlator to the exponential time dependence. 

In the case of the flavor singlet correlator, the lightest state in the sum (jS)) is the vacuum 
state. Its corresponding coefficient cq (jH) is the scalar condensate (qq). Another important 
light state that contributes at large t is vrvr, so extraction of rricr requires the fit to 

c{t) * Id + YlhX^-''^^^' + i^^)" ■ (10) 

The extraction of m„ is very challenging since C{t) requires the calculation of the discon- 
nected diagram (see previous Section) and since RHS in (fTUI) is largely dominated by (qq)^. 

These two problems do not affect the study of the flavor non-singlet meson. However, 
even in this case there are several multi-hadron states which are light and need to be taken 
into account in the fit of the correlator Q at large t in order to extract mg. The lightest 
multi-hadrons states with = O"*" are two-pseudoscalar states in S'-wave. In case of / = 1 
correlator, the contribution of scalar meson oq is accompanied by contributions of vrr], KK 
and ITT]' in three-fiavor QCD. Let us note that in nature these three states are lighter than 
observed resonance ao(1450); the state irrj is also lighter than observed resonance ao(980). In 
two-fiavor QCD, the only two-pseudoscalar state nt]' is relatively heavy and not so disturbing 
for the extraction of niao from ©. 

The above derivation of time-dependence for a correlator was based on QCD, which is a 
proper unitary field theory. The resulting correlator © is positive definite. Let us point out 
that certain approximations used in lattice simulations (quenching, partial quenching, stag- 
gered fermions, mixed-quark actions) break unitarity and may render negative correlation 
function. These approximations will be discussed in Section 5 together with the necessary 
modifications of the fitting formula Q. 
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4 Mass of scalar meson with 1=1 



A lattice simulation of the scalar meson oq with J = 1 ^ is presented in this section, as an 
example. It employs two dynamical quarks^, lattice spacing 0.12 fm, lattice volume 16^ x 32 
and ensemble of about 100 gauge configurations jH |3] . The advantage of simulation j3] is 
that its discretized (Domain- Wall) fermion action has good chiral properties: it is invariant 
under the chiral transformation for rriq = even at finite lattice spacing'^, which is not the 
case for some of the commonly used discretized fermion actions. Another advantage of the 
simulation with two dynamical quarks |3] is that the exponential fit of the correlator at 
large t renders rriao- The conventional exponential fit is justified in this case since the only 
two-pseudoscalar intermediate state in two-flavor QCD is nt]', which is relatively heavy and 
does not affect the extraction of (see previous Section). 

The resulting mass is presented in Fig. |21for different input masses ruu^d, where isospin 
limit rriu = nid is employed. There are no simulations at physical masses rriu^ since the pion 
cloud around the scalar meson with = /ic/140 MeV~ 9 fm would be to squeezed on the 
lattice with extent 16 x 0.12 fm~ 2 fm. The u/d quarks and pions are heavier in simulation 
than in the nature in order to avoid large finite volume effects. The linear extrapolation of 
to the physical quark mass ntu^d — 4 MeV in Fig. IHl gives 

m^o = 1.58 ± 0.34 GeV . (11) 

Although our result for the mass of the lightest qq state with 1=1 has sizable error-bar, it 
appears to be closer to the observed resonance ao(1450) than to ao(980). It gives preference 
to the interpretation that ao(980) is not conventional qq state. 

Results from other lattice simulations of the light scalar mesons can be found in [B]-|llj. 
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Figure 3: The triangles present resulting mao for three values of bare quark masses niu^d HI- The 
dashed line is the linear extrapolation of mao to the value of mu,d in nature. 



^Fermion determinant in incorporates quarks i = u,d. 

^This is strictly true only when the 5th dimension in Domain- Wall fermion action is infinitely large. 
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5 Problems due to unphysical approximations 



The simulation presented in the previous section is a discretized version of two-flavor QCD 
and does not employ any unphysical approximations except for the discretization of space- 
time. It renders positive definite correlation function, as expected in proper Quantum Field 
Theory ©. 

However, lattice simulations often employ unphysical approximations which facilitate 
numerical evaluation. One of the indications that the simulation does not correspond to a 
proper QCD is the negative scalar correlator. Another sign of unphysical simulation is when 
1 = 1 correlator drops as e~^*^'"* at large t although the lightest two-pseudoscalar state with 
/ = 1 is nr]. Both of these unphysical lattice results can occur if the theory that is being 
simulated is not unitary, which is the case for all the commonly used approximations listed 
below: 

• In gwenc/ied simulation the fermion determinant in Q is replaced by a constant. This 
corresponds to neglecting all the closed sea-quark loops. The 1 = 1 scalar correlator 
is negative in this case and its negativity was attributed to the intermediate state 
TTT]' in Ref. [Oj. The prediction for irrj' intermediate state in quenched version of 
Chiral Perturbation Theory (ChPT) describes the sign and the magnitude of the lattice 
correlator at large t well jHl Ej. The mass m^o was extracted [HI EI by fitting the 
quenched 1 = 1 correlator to the sum of e""^""* term and the contribution of irt]' as 
predicted by Quenched ChPT. 

• In partially quenched simulation the mass of the sea quark is different from the mass of 
the valence quark, although they are the same in nature. The mass of the valence quark 
is the mass that appears in the propagator of the correlator Co (0), while the mass of 
the sea quark is the mass that appears in the fermion determinant (jH). The partially 
quenched scalar correlator with 1 = 1 was found to be negative if rrivai < fT^'sea This 
was attributed to intermediate states with two pseudoscalar mesons and was described 
well using partially quenched version of ChPT The mass rriaQ was extracted by 
fitting the partially quenched correlator to the sum of 6"™""°* term and the contribution 
of two-pseudoscalar states as predicted by Partially Quenched ChPT jlj. The resulting 
mass agrees with the mass (fTT|) . 

• The simulations with mixed quark actions employ different discretizations of the Dirac 
operator for valence and sea quarks. The method of extracting scalar meson mass from 
a such simulations was proposed in [T^ IT^ . 

• The simulations with staggered quarks use an artificial taste degree of freedom for 
quarks in order to solve fermion doubling problem [3]. The method of extracting 
scalar meson mass from simulations with staggered quarks jJJ was proposed in [T^ . 

All these approximations modify the contribution of two-pseudoscalar intermediate states 
with respect to QCD. The effects of these approximations can be therefore determined by 
predicting the two-pseudoscalar contributions using appropriate versions of ChPT. These 
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analytic predictions [HI IH El El allow the extraction of the scalar meson mass from the 
correlator as long as the contribution of two-pseudoscalar intermediate states does not com- 
pletely dominate over the e~™^* term. 

6 Conclusions 

The nature of scalar resonances below 1 GeV is not established yet. A lattice determination 
of the masses for ground qq scalar states would help to resolve the problem. 

In principle, the scalar mass can be extracted from the scalar correlator that is com- 
puted on the lattice. However, the interesting term e""''^* in the correlator is accompanied 
by the contribution of two-pseudoscalar states e~^^^*. The problem is that the energy 
of two-pseudoscalar states is small, so they may dominate the correlator and complicate 
the extraction of scalar meson mass. On top of that, the contribution of two-pseudoscalar 
states is significantly affected by the unphysical approximations that are often used in lattice 
simulations. Luckily, these effects can be predicted using appropriate versions of Chiral Per- 
turbation Theory and they agree with the observed effects on the lattice correlators. We give 
the list of references, which provide the expressions for extracting ms from the correlators 
for various types of simulations. 

A simulation, which does not suffer from the problems listed above, gives 1.58 ±0.34 GeV 
for the mass of the lightest qq state with 1 = 1. This supports the interpretation that 
observed a0(1450) is the lightest {qq)i=i state, while a0(980) might be something more 
exotic. 
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